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Abstract 

We propose a low energy effective theory of lightly doped two-leg t-J lad- 
ders with the help of slave fermion technique. The continuum limit of this 
model consists of two kinds of Dirac fermions which are coupled to the 0(3) 
non-linear sigma model in terms of the gauge coupling with opposite sign of 
"charges" . In addition to the gauge interaction, there is another kind of attrac- 
tive force between these Dirac fermions, which arises from the short-ranged 
antiferromagnetic order. We show that the latter is essential to determine 
the low energy properties of lightly doped two-leg t-J ladders. The effective 
Hamiltonian we obtain is a bosonic Gaussian model and the boson field ba- 
sically describes the particle density fluctuation. We also find two types of 
gapped spin excitations. Finally, we discuss the possible instabilities: charge 
density wave (CDW) and singlet superconductivity (SC). We find that the SC 
instability dominates in our approximation. Our results indicate that lightly 

doped ladders fall into the universality class of Luther-Emery model. 
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I. INTRODUCTION 



The properties of ladder systems described by t-J and Hubbard models have been the 
subject of intensive studies recently PHTB(]. The reason is that there are systems such as 
(VO) 2 -P207 0] and SrCu 2 3 JT^] which can be the possible realization of these models. 
Experiments on magnetic susceptibility and neutron inelastic scattering show the existence 
of a finite spin gap. Moreover, a recent measurement shows the sign of superconductivity 
in the doped spin ladders Therefore, the study of these systems can offer new insights 



into the nature of magnetic, charge density, and pairing correlations in strongly correlated 
electron systems. 

The undoped ladder systems show unusual magnetic behaviors. When the number of 
chains is even, there is a spin gap. When it is odd, the spin excitation is gapless. This can be 
easily understood by considering the limit of strong rung interactions. For two-leg ladders, 
since the spins on every rung forms a singlet first in this limit, it can be considered as a set 
of weakly coupled rung singlets. The spin excitation is formed by turning over one spin on a 
rung and this costs finite energy. Along the same reasoning, a three-leg ladder is effectively 
equivalent to a spin-1/2 chain and the latter has gapless spin excitations. Numerical studies 
|17|| show that the spin gap for two-leg ladders persists even at the experimentally interested 



isotropic point. (By this we mean that the rung interaction is equal to the intrachain 
interaction.) This indicates that even at the isotropic point, the rough picture of a ground 
state dominated by rung singlets is robust. People have investigated spin ladders by the 
semiclassical (large spin) approach |nj. Remarkably, this approach qualitatively captures 
the basic feature of this system as in the case of spin chains. Besides, the spectrum predicted 
by the non-linear sigma model is in good agreement with other approximations |19[ both 
in the strong and weak rung interaction limit. These facts imply that the non-linear sigma 
model correctly describes the low energy sector of the spin liquid state of ladder systems 
though it ignores fluctuations along the rung. This is one of the motivations that we would 
like to use the large spin approach to study the lightly doped case. 

The effect of doping on the antiferromagnetism is an important and unsettled issue 
in strongly correlated electron systems. The main difficulty lies in hole motions, which 
cause frustration of the antiferromagnetic (AF) order. Theoretically, we have no adequate 
analytical methods to deal with the competition between charge and magnetic fluctuations 
because there is no obvious small parameter which facilitates an expansion about a tractable 
model. In the present paper, we employ the large spin expansion which permits us to map 
the original model to a continuum field theory. And this allows us to tackle the problem 
with analytical methods due to the 1-d nature of this system. Mean field studies [fj-f|] 
and numerical calculations show that the spin gap in undoped two-leg ladders still 

persists at low doping concentration. This implies that the underlying short-ranged AF 
order is not destroyed too much by hole motions as the doping concentration is low. Thus, 
after we resolve the constraint in the slave fermion representation, we assume that spin 
variables have a strong short-ranged AF order. This leads to the t — J model proposed 
by Wiegman, Wen, Lee, and Shankar p0|,2l|. Because of this background AF order, it 



is natural to say that there are two kinds of holes (on A and B sublattice) with opposite 
sign of (fictitious) charges and they are coupled to the staggered magnetization through the 
Berry phase. In contrast to previous studies [PU |, we keep the nearest-neighbor attractive 
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four-fermion interactions between A and B holes. This attractive force can be understood 
as the following: the energy of two holes on the same spin singlet is lower than that of 
two holes on different singlets because there are less broken bonds in the former. And 
this is equivalent to an attractive force between holes on different sublattices. We discuss 
this model in the absence of the quartic fermion interaction first. To study the low energy 
physics we can linearize the dispersion relation of fermions about the fermi points . It turns 
out that we have four branches of massless Dirac fermions coupled to the non-linear sigma 
model. From the study of Schwinger model, we know that massless Dirac fermions will 
screen the long range Coulomb force or give a mass term to the gauge field on account of 
the chiral anomaly. Therefore, the gauge fields are in the Higgs phase ||22|| . Because of 
the Gauss's law, all excitations must be gauge singlets (or mesons in the jargon of particle 
physics). In the present case, we have three gapless spin-singlet collective modes. There are 
also gapped spin-1/2 excitations which carry the electronic charge but no magnons. The 
SC instability is enhanced. Then we switch on the nearest-neighbor attractive four-fermion 
interactions. Two of the gapless modes become massive and the long range Coulomb force 
is not screened, i.e. the gauge fields are in the confining phase. Consequently, we have only 
one gapless charge mode, which is spin-singlet. We also have gapped spin-triplet excitations 
and electron-like collective modes, which are spin-1/2 and carry electronic charges. We 
compute the exponents of pairing correlation function and Akp CDW susceptibility. In our 
approximation, turning on the quartic fermion interactions enhances the tendency toward 
SC relative to the tendency toward CDW. Therefore, a weak interladder interaction will 
lead to SC, which has been predicted by other approaches 0-f|,[7|-|9|] . Our analysis indicates 
that these two types of attractive forces play different roles in two-leg t-J ladders. While 
the formation of spin gap and spin-hole bound states are due to the gauge interaction, four- 
fermion interactions between A and B holes are responsible for the pairing between holes. 
In addition, inclusion of the latter drastically changes the low energy properties of doped 



ladders such that they fall into the universality class of Luther- Emery model p3fl . Our 
results confirm the conclusions from numerical studies |J. 

The rest of the paper is organized as follows. In section two we derive the low energy 
effective action. In section three we discuss the implications of this action. We study the 
effect of gauge interaction in subsection A. In subsection B we take into account the four- 
fermion interaction between A and B holes. Section four is the conclusion. We give a 
summary of the bosonization rules we used in the appendix. 



II. DERIVATION OF THE EFFECTIVE ACTION 



Since the first part of our derivation is valid for general t-J models, we will not write 
down the ladder index explicitly until it is necessary. We start from the following model 

H = -t ^X a0 (i)X 0a (j) + $>X ff(7 (z) + 

{i>j) 
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where means the nearest neighbor sites, X ab =| a)(b | and | a) =\ 0), |f), ||) corre- 

sponding to the empty site and spin-up (spin-down) sites. Since transitions between empty 
and occupied states include a change in the fermonic number, the operators X a0 (i), X Qa (i) 
are fermonic and the operators X a(T >(i), X 00 (i) are bosonic. It is easy to check that they 
satisfy the following graded Lie algebra called Spl (1,2) ]23J: 

[X ab (t),X cd (j)}± = 5 t] (5 bc X ad (i) ± 5 ad X ch (i)) (2) 

where (+) should be used only if both operators are fermonic. In the limit of zero doping, 
e d — > — oo and the model is reduced to the spin-1/2 Heisenberg model. 

Among various methods to deal with eqs.(|I|) and (Q), there are two popular representa- 
tions of the above graded Lie algebra - - slave fermion and slave boson. We adopt the former 
and introduce a vacuum state annihilated by operators fi and b a (i). Then the X-operators 
can be represented as follows 

X Q<T (i) = f?K(i), 
X a0 (i) = f i b+{i), 

x aa ,(i) = K(i)K>(i), 

Xoo® = ftfi (3) 

with the constraint b^(i)b a (i) + ftfi = 1 and o = 1, —1 for spin up and down, respectively. 
Here fi, ff satisfy canonical anticommutation relations and b a (i), b^(i) satisfy canonical 
commutation relations. To use the large spin expansion, we replace the above constraint 
with this one: b+(i)b a (i) + ft fi = 2S, which is called the spin-s representation of Spl(l, 2). 
Now the X-operators represent transitions between states | s) and | s — 1/2). With the help 
of eq.(H), we can write down the path integral representation of the t — J model 

Z = J D[f]D[f + ]D[b a ]D[b + a ]5(K^(i) + ft ft - 2S)exp(-I), 

1= f dT^iftidr- M + K(i)drh(i)]+ f dTH, 

Jo i Jo 

H = tJ2ftfri(*)bAj) + 

{ J2[bt(i)h(i)bi(j)b a (j) - lbt(i)b a (i)b+(j)bfi(j)}- (4) 

In the above eq., e d has been absorbed into /i, the chemical potential of holes in terms of 
the constraint in the path integral measure. Based on the property of Grassmann variables: 
(f + f) 2 = 0, we can resolve the constraint in the measure |24[] as the following 




where S = (S — |/ + /)fi and fi = (sin 9 cos (j), sin 6 sin <fi, cos 6) . With the new variable ft, 
we can rewrite the partition function as follows 
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Z = J D[n]D[f]D[f + ]exp(-I), 



j 

r* 3 

+ / dr H, 
Jo 

j Ea - ^tf/oa - ^./y • «i (6) 

where t = 25* and 7 = 4JS* 2 . (Note that when S = 1/2, t = t and J = J.) (fii | fi 2 )s 
represents the overlap of two spin-S" coherent states, | Qi) and | S~2 2 ) • A(f2) is the monopole 
vector potential, which satisfies y x A = ft. (Here we have chosen the gauge such that 
A ■ d T fl = — cos 9d T (j) and x is independent of r.) 

Now we focus ourselves on two-leg t — J ladders. Using the notation of eq.(|5]), the action 
can be written as follows 



f dr £ £[#m(3r - ,/)/,,„. + i(S - - y fi,J h , 

J ° m j z 

f/3 

A(n j>m ) ■ d T ft jtm \ + J dr H, 



H = t Yjftmkmi^hm | %m)l/2 + *_L YX&oA^J,! I n i,2>l/2 + h.C.) 
m (ij) 3 
J 1 1 

+ J -f EC 1 - - • «^ (7) 

where m — 1,2 is the chain index. To proceed, we have to make some assumptions. First 
of all, we assume that there is a short-ranged AF order such that we can parameterize Jjm 
as the following |Tj| 



n i>m = (-iy +m y/l - a 2 L^n(x) + aL m (x) (8) 

where n 2 = 1, n • L m = and a is the lattice spacing, n is the order parameter and L m are 
the fast modes which will be integrated out. Comparing the £-term and J-term in eq.(^), 
it is clear that in the large spin limit J-term dominates and thus eq.fl8|) is valid. For small 
spins, we expect it is still a good assumption if there are strong short range AF correlations 
and the latter is reflected by the existence of a finite spin gap in two-leg ladders. However, 
some conditions are required for a non- vanishing spin gap. Firstly, according to numerical 
studies 0-0, the hole concentration must be low. Secondly, J/t should be of order one. 
As has been noticed in ref. 0, the Nogaoka theorem is applicable in ladder systems. This 
theorem says that the one-hole ground state is ferromagnetic at J = 0. This phase may 
be stable when J/t is very small fill . Besides, numerical studies show that the phase 
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separation occurs when J/t > 2. Therefore, the qualitative feature of the following results 
may be valid for the spin one-half case when the ratio of J/ 1 is order of one and the hole 
concentration is low. 

Because of the background AF order, the original lattice is divided into two sublattices 
and it is natural to have two kinds of holes. We will call them A-holes and B-holes when j+m 
is even and odd, respectively. An immediate consequence of eq.(§) is that the hole cannot 
hop coherently between different sublattices, i.e. the intersublattice hopping is forbidden. 
This can be understood as the following: The amplitude for a hole hopping from site i to 
j is the product of the overlap of spin and orbital wavef unctions, i.e. — ty (fij | fij). In this 
case, it is —t(flA \ £Ib)- In view of eq.(^), this is zero if we have perfect Neel order (recall 
that (f2 | — f2) = 0) and exponentially small in a system with strong short-range AF order. 
Thus, t- and tj_-terms are effectively removed from the low energy effective Hamiltonian in 
the large S analysis. One of their effects is to renormalize the parameters in low energy 
physics. It seems that the hole can hop coherently within the same sublattices under the 
hypothesis eq.(f§). However, the question of whether the hole can have coherent hopping is 
related to whether we have well-defined quasiparticles in low energy. As pointed out in ref. 



25| , the latter depends on the density of states of spin waves in low energy. It turns out 
that there can be no well-defined quasiparticles in one dimension if the spin excitation is 
gapless. Fortunately, there is a spin gap in our system. Following the calculations in ref. 



|25| , there will be a sharp peak at the energy scale —t in the spectral function of one-particle 



Green function and this implies a well-defined quasiparticle band at the bottom of the hole 
spectrum. Consequently, via the interactions with spins, the holes are able to acquire a 
kinetic energy, whether there is a "bare hopping term" or not. In the following, we shall 
introduce a t'-term, i.e. hopping within the same sublattice, to represent the kinetic energy 
of holes. According to numerical calculations , t is of the same order as J. 

Substituting eq.(§D into eq.(^), adding a gauge- invariant t'-term to eq.(^) and expanding 
the action to the quadratic terms of L m then integrating out L m , we get to the following 
effective action 

I = I n + I h , 

I rP r 1 
/„ = — ^ / dr dx (— | d T n \ 2 + | d x n | 2 ), 



2g 2 Jo ' " ,2 



rp 

J m j 

ftMKdr + ia - »)fB,mti)} + f dr H h (9) 

Jo 



and 



Hh = H + Hi + H 2 , 

H =^ E HifimU + 2)fA,mU) ^P KAx) + h.C.] + 
m j 

(A -> B,a x -> -a x ), 

H l = -J EEi : fXm(j)fAAj) ■■■ ftnti + l)fB,mU + 1) = +(^ - B)], 
rn j 



6 



H2 = —f Et : />.;!./)/>.:(./) " f^fsM : + (A - B)\ (10) 

i 



where # 2 = 2/Ja(l - o"/ 2 ^) 2 , w s = g(l - 6/2S)Jl + fy, t' = 2St' , Ax = 2a, and 5 is the 
hole concentration. In the derivation of above eqs., we have used a property of spin coherent 
states: (f2i | 2 )s ~ exp — iSA ■ (Hi — f2 2 ) when f2i ps f2 2 - The gauge fields = |A • <9 M n. 
It is clear that A and B holes carry opposite sign of charges because they are on different 
sublattices in which the staggered magnetizations have opposite directions. We propose that 
the low energy sector of lightly doped two-leg ladders can be described by eqs. (Q) and (fLOD . 

We conclude this section with a brief summary of the CP 1 representation of the non-linear 
sigma model |26fl , which will be used later. The order parameter n can be parameterized 
with a normalized spinor as 

n = zaz , zz — 1. 

The CP 1 version of the non- linear sigma model in the Euclidean space is 
Z = J D[a^]D[z\D[z]D[\}exj){-I), 

I = i J d 2 x [| 8^z - ia^z | 2 +iA(| z | 2 -1)] (11) 

where xq = v s t, e 2 = 2g 2 v s and A is the Lagrangian multiplier. An appropriate choice of 
the gauge makes a M in eq.flTTD equivalent to the ones in eqs.fl^) and ([TDD an d thus we can 
use the same notation to denote them. In the CP N model at large N, the ^-quanta become 
massive, i.e. A acquires a non-vanishing mean value. At distances larger than this scale, 
the effective action of gauge fields will have the usual Maxwell term. Accordingly, the low 
energy excitation is a massive triplet, which can be considered as the bound state of z and 
z. 



III. EXCITATION SPECTRUM 



To discuss the low energy physics, we linearize the dispersion relation of fermions about 
their fermi points, which satisfy kpa = |(1 — 8). ( We have assumed that the Hamiltonian 
is invariant when we interchange the chain index and A and B holes. Besides, remember 
that a is the lattice spacing of the original lattice.) After doing that, Ho becomes 

H = E + v Y^ / dx [: ^ m a{-id x - a x )ifj A ,m ■ +(A -> B, a x -> -a x )\ (12) 

m J 

where vq = At' asm.2kpa = 4t'asin7ro" is the fermi velocity, ip is the two-component Dirac 
fermion, and a = a^. Now we can analyze the effects of different interactions. 
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A. The Effect of Gauge Coupling 



We set Hi = H 2 = first. Then we rescale the imaginary time: Vqt — > r and do 
analytical continuation to the real time formalism. The effective action is 

I = I + In, 

= £ / d2x U>A, m Y{idu eA^ipA^ + 

(A-> B,a u ^ -a M )] (13) 

and J n has the same form as the one in eq.(Q) and 7 M is the Dirac 7-matrices. Here A n 
are the external electromagnetic fields. In terms of the standard bosonization rules (see the 
appendix), I can be bosonized as the following 



J = \ £ / # X l(d^A,m) 2 + {d^B, 

"^=£ / A^ V d v {(j) A>m + (j)B, m ) + 

y/TV m J 

A=E I SX a ^ U dv{^m ~ (f>B,m). (14) 
V 71 " m J 

We define 0±, m = ^(0A,m ± 4>B,m)- Then J can be written as 

^ m J 

\/f e E/ ^ V"^+, m + yl £ / ^ a^d v <l>-, m . (15) 
The above action can be further simplified if we define the following fields 



$1 = — (0+.1 + 0+.2), 



1 

$2 = - 0+,2), 



*4 = -^(0-,l-0-,2). (16) 

With the help of the above canonical transformation, we obtain the following low energy 
effective action 



*> = -£/ ^ (^$ a ) 2 - -J=e J d 2 x A^d u ^ 

J d 2 x ^ v d^a v . (17) 



a=l ' 
2 



What can we learn from eqs.fllTD and (|T7|)? Firstly, there are three gapless spin-singlet 
excitations: $1, $2 and $4. Secondly, to understand the coupled system of $3 and z, we have 
to investigate the gauge field dynamics. This can be done by integrating out the $ 3 -field 
since its action is quadratic and we obtain 



— J d 2 xd 2 y f(x)ln \x-y\ f(y) 



where / = e^ u d^a u is the dual field strength. If we choose the Lorentz gauge: = 0, 

then the above eq. becomes the mass term of a^. That is to say, the fluctuations of the 
$ 3 -field will screen the long range Coulomb force. It follows that there are massive spin-1/2 
excitations and they are neutral with respect to the a^-fields in view of the gauge field mass. 
(This is easy to be seen by integrating Gauss's law from —00 to +00.) Since the ^-quantum 
has charge 1, where is the compensating charge from? It must come from the hole sector. 



We can understand this as the following pi 



Let us integrate out all fields except $3. The action we obtain must be of the form 
cosS-y/vrn^ where n is an integer. (Of course, the renormalization of the kinetic term is 
possible.) This is because the coefficient of $3 in eq.(|l7|) is proportional to the instanton 
density: -r-e^ v d^a v . The action is invariant under the transformation: $3 — > $3 + y/n/4: 
and the cosine just has the form which satisfies the requirement. If the cosine is relevant, 
this translation symmetry is spontaneously broken and the value of $3 is pinned at some 
minimum of the potential. In addition, its solitons carry exactly the charge needed to 
match the charge of the z-quantum. The Gauss's law demands that at least one of the 
cosines in question must be relevant and thus the excitations corresponding to the $3 sector 
are massive. These solitons are nothing but the A (or B) holes. They carry the electronic 
charges as well as the U (1) charges. As a consequence, we expect these massive doublets also 
carry the electronic charges, i.e. they carry the same quantum numbers as electrons. We 
should emphasize that there are no zz bound states, i.e. massive triplets, in this model. As 
we shall see later, the absence or presence of this type of excitations is one of the distinctions 
between the Higgs phase and confining phase in our model. 

After identifying the spectrum, we would like to examine the pairing correlation function. 
Because of the spin gap, the pair-field of singlet SC can be defined as the following 

A(j) = f A ,iU)fBM- (18) 

Since the correlation functions of $3 decay exponentially, the long distance behavior of the 
pairing correlation function can be calculated by the following operators: 

A(j) ~ exp (-iy/irQx) exp (-iy/TT& 2 ) exp (-2^6)4) + 

exp (— i^/nQi) exp {1^/1:^2) exp (— 21/^64) (19) 

where G a is the dual field of $ a . The pairing correlation function behaves as 
(A(j)A+(0)> j-^. Compared with the one of free fermions, which behaves like pW, 

we see that the pairing susceptibility is enhanced. 

To sum up, if we consider the gauge coupling only, then the low energy effective Hamil- 
tonian consists of three gapless spin-singlet modes and one massive spin-1/2 modes, which 
carry the electronic charges. (It is possible that there are massive spin-singlet modes. How- 
ever, the gapless modes and the massive doublet predominant the long distance behavior of 
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correlation functions.) Compared with the numerical results 0, there are too many gap less 
modes but no magnons. It is not enough to consider the gauge interactions only. We will 
see in the next section the importance of taking into account H% and H 2 to get the correct 
low energy properties. 



B. The Role of H ± and H 2 

Here we take Hi and H 2 into account. The continum limit of them are as follows 

m J 

g X COS TX5 ^ dx [lpA,L,m' l pA,R,m1pB,R,m' l l J B,L,m + (R L)], 
m J 

H 2 = -92 J dx [: VXiV^M " ^b^b.i ■ +(A <-> B)] - 

g 2 J dx [^A,R,l^A,L,l^t,L,2^B,R,2 + ^A,L,1^A,R,1^,R, 2 ^B,L,2 

+(A <-> B)} (20) 

where g\ = 2Ja/v , g 2 = J±a/v and ipL,R are left-handed and right-handed fermions, 
respectively. With the same convention used in the previous section, the bosonized forms of 
the above eqs. are as follows 

h x = -I 1 E / dx l(d x <i> + , m ) 2 - (o>,0-, m ) 2 ] + 

9\ / dx cos \/87T0_ jm , 
H 2 = -— dx (d x (f) +)1 d x (f) +y2 - d x (f> id x <f>-^) - 

71 J 

2g 2 J dx cos v / 2^r(0 +i i — </) +i 2) cos v / 27r(0_ i i + 0__2). (21) 



We can do the following canonical transformation to diagonalize the quadratic part of the 
Hamiltonian 

$1 = \/^(0+,l + 0+,2), 

$2 = y^(0+,i - ^+,2)> 

$3 = ^(0-,l + 0-,2), 

$4 = ^(0^1-^,2). (22) 
The parameters in the above eqs. are as follows 
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1 


9i + 92 


71 


1 


9i - 92 


7T 


1 + 


9i + 92 


7T 


1 + 


9i ~ 92 


71 



(23) 

These relations are valid only in the weak coupling limit. In our case, this corresponds to the 
large spin limit. (Remember that t'(= 2St ) appears in the denominator of the definition of 
gi.) We can see that K\ < 1 and K 3 > 1 because <fc > 0. (This corresponds to the attractive 
force between holes.) The effective Hamiltonian of the hole sector is 



Hh = \ I>a / dx i(d x Q a ) 2 + (<9 X .$ Q ) 2 ] + 

a=l 

2g 1 / dx cos J — $ 3 cos J — <J> 4 - 
j V As V Ka 



i f / 47T / 47T 

2# 2 / da; cos W— $2 cos W— $ 3 - (24) 

J V A 2 V -^3 

Here, we neglect the gauge coupling temporarily and will discuss it later. The relevancy of 
these interactions is determined by their scaling dimensions. They are = + and 

A 2 = + which correspond to and g' 2 -term, respectively. From eq.(p3|), K 2 < 1 
and > 1 at the isotropic point. This implies Ai < 2. Therefore, g^-term is a relevant 
operator in the sense of renormalization group (RG). $3- and $4-field are pinned at some 
values and both acquire gaps. They are decoupled from the low energy theory. Taking into 
account these facts, the effective Hamiltonian becomes 



H h = - K * J dx [(d x Q a ) 2 + (d x $ a ) 2 ] -g J dx cos J— $ 2 (25) 

J a=l,2 " V 2 



where g = 2g 2 (cos yj^$3)- In the large spin limit, K 2 > 1/2. Thus the scaling dimension 
of g-term is less than two. It is a relevant operator. The $2-field is also massive. In the low 
energy limit, we obtain our effective Hamiltonian as the following 

He ff = ^- J dx[{d x Q x f + {d x ^ x f]. (26) 

Now we would like to discuss the implications of our results. First of all, the low energy 
effective Hamiltonian consists of a gapless spin-singlet mode and from eq. (|22"[) , this mode 
describes the total charge density fluctuation. Furthermore, there is only one free parameter 
K\, the compactification radius of $i-field, which has to be determined from experiments. 
This supports the suggestion proposed in ref. f|. Second, since the <3>3-field is pinned at 
some value, it can not affect the non-linear sigma model too much. Especially, it is unable to 
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screen the long range Coulomb force. As a result, the gauge field is in the confining phase. 
This leads to two types of spin excitations. One is the massive triplet (In the CP 1 language, 
it is the zz bound state), which is the same as the magnon in the undoped case. The other 
is the bound state of holes and z quanta, which carries the same quantum numbers as the 
electron. In fact, this excitation can be considered as the breaking of a hole pair and we 
have two quasiparticles. Each carries electric charge one and spin one-half. The latter is 
also observed in ref. ||. In that paper, Troyer et al. find that the spin gap is determined 
by the bound state of spinons and holons. We cannot compare the gaps of these excitations 
in our approach. However, we give a picture about the formation of the unexpected spin 
excitation. The Berry phase term gives rise to the necessary attractive force between spinons 
and holons. Since this force is a gauge interaction and the latter in Id is a linear confining 
potential, this results in bound states of spinons and holons. 

After we understand the spectrum, we can calculate the asymptotic behavior of various 
correlation functions. The most important ones are pairing and CDW correlation function. 



The definition of the pair- field is the same as eq.([18|). The long distance behavior of the 



pairing correlation function can be calculated by the operator 



A(j)~exp(-z v /7rir 1 e 1 ) (27) 
because other fields are massive. The result is 



(A(j)A+(0)) ^ (28) 



Next we shall examine the 2k p CDW susceptibility. The corresponding order parameter 
Ocdw can be expressed as the product of hole operators and spin operators. As emphasized 
in the previous paragraph, matter fields do nothing much on the spin sector of the effective 
action. We can find the contributions of spin operators from results of the undoped case. 
The work of Shelton et al. |?7| showed that the low energy theory of two-leg ladders is 
described by four decoupled non-critical Ising models with three of them in the ordered 
phase and one in the disordered phase (or vice versa). Moreover, the spin part in Ocdw can 
be expressed in terms of the order (a) and disorder (/i) parameter fields of the Ising model as 
^1^2/^3/^0 or <7i<72cr 3 cro- It is clear that the correlation functions of the above operators decay 
exponentially to zero. This implies that (Ocdw(x)Ocdw(®)) shows the same behavior. 
Then we have to consider the 4fcp CDW susceptibility or (O CDW (x)O CDW (0)). Using the 
OPE ii(z)fj,(u) ~ l/(z — uj) 1 I a + • ■ -, it is straightforward to see that the spin part of this 
correlation function contributes a non- vanishing constant to it. Thus, the long distance 
behavior of (O CDW (x)O CDW (0)) can be determined solely through its hole part. If we 
define the operator 



Oh(x) = TpRAl^B^LA,^ 



L,B,2 



then 



/ 4ir 

exp(- M /— $0, (29) 

V K i 



(O 2 CDW (x)O CDW (0)} ~ (O h (x)O h (0)) 

~ uk- (30) 
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In the second line of eg. (p9|) , we keep the gapless mode only. From eqs. (|28|) and (p0|), it 
is clear that SC dominates when K\ < 2 and CDW dominates when K\ > 2. In the 
large spin limit, Ki « 1. Therefore, we conclude that SC susceptibility dominates in two- 
leg t — J ladders and a weak interladder interaction will lead to superconductivity at low 
temperature. Also the exponents of pairing and CDW susceptibility satisfy the relation: 
Ksc ■ Kcdw = 1- We arrive at the same conclusions as previous numerical investigations. 
Although we heavily rely on the large spin approximation, our results should capture the 
basic feature of the system with spin one-half. 

In summary, after we take Hi and if 2 into account, the low energy effective Hamiltonian 
only consists of one gapless charge mode, which describes the charge density fluctuation. The 
spin excitations are electron-like quasiparticles and magnons and both have energy gaps. The 
2/ci? CDW susceptibility shows exponentially decaying behavior while those of 4kp CDW and 
singlet SC both show power-law behavior. With the above results, we conclude that lightly 
doped two-leg ladders fall into the universality class of Luther-Emery model. In addition, 
this phase is dominated by singlet SC susceptibility according to our analysis. 

There are two related works which should be mentioned. The first is the paper by 
Ichinose and Matsui [O]. They also adopted the slave-fermion and CP 1 boson technique 



to treat the two-leg t — J ladder. The other is the work of Ivanov and Lee [|13j. They use 



the slave-boson scheme and arrive at the same conclusion that the low energy sector of the 
two-leg t — J ladder is a Luther-Emery liquid. For quantitative description of this system, 
they introduce a new order parameter n pa i r = (rii — 5) (ri2 — 5) where rij is the hole density in 
the i-th chain in addition to the ones we considered. They suggested that the long distance 
behavior of the 2kp part of (n pair (x)n pair (0)) (Notice that their 2kp is equivalent to our 
4k p.) is as the following 

A B 
+ 



X a ^ 

where A and B depend on the average pair overlap: A ^> B at £ pa i r 5 <C 1 and A < B at 
ipair& 1- Here £, pair is the size of the hole pair and is of the order of the lattice spacing. 
a± is the exponent of (Oq DW (x)Oq DW (0)) . We calculate this correlation function and we 
obtain a power-law behavior with the exponent equal to 2/Ki, which is the same as the one 
of (Oq DW (x)Oq DW (0)) . We do not find an exponent equal to 2/i£i+2 even including higher 
harmonics. Since our approach is valid when 5 <C 1, we cannot tell whether the second term 
exists or not at larger doping concentration. 



IV. CONCLUSION 



Anderson has proposed that the spin liquid state may evolve into a superconductor upon 
doping. However, it has been proved notoriously difficult to have any concrete analytic 
result to confirm this idea in 2D. A doped spin ladder may provide a good place to study 
this mechanism though there is no true long range order here. This is so because the undoped 
two-leg ladder is a kind of spin liquid state and it is simpler to deal with this problem both 
analytically and numerically. Soon after the discovery of high T c superconductivity a model 
was proposed to describe the doped spin liquid state by coupling holes to the non-linear sigma 



13 



model and the gauge interaction between A and B holes provides the necessary attraction 
for pairing. Our approach is basically to apply the same idea to ladder systems. In contrast 
to previous studies, the gauge interaction plays a minor role on the formation of hole pairing 
in two-leg t — J ladders. Instead, it is mainly due to the effective attraction between the 
nearest-neighbor holes. Therefore, upon doping, these holes are inclined to stay on the 
same preexisting singlets. Although this may be a unique characteristic of two-leg ladders 
exclusively, this point deserves further examination in higher dimensions. We also point out 
that the gauge interaction is responsible for the electron-like collective mode found in ref . || . 
The existence of this type of excitations is independent of the quartic fermion interactions. 
The absence or presence of the latter determines whether the long range Coulomb force is 
screened or not and thus affects the magnitude of the gap of excitations with the non-trivial 
spin quantum number. More importantly, it does affect the existence of the magnon: there 
are no magnons without the quartic fermion interactions. We have to emphasize that both 
the effective attraction between nearest-neighbor holes and gauge interactions are due to the 
strong short-range AF background, which has its origin in the strong repulsive interactions 
between electrons. 

The main assumption we made is to replace eq.(|7|) with eqs.([|) and (|i~0p. This is justified 
for undoped case and for lightly doped two-leg ladders. The spin excitations in the undoped 
three-leg ladder are gapless and thus the distortion of the short-range AF background arising 
from hole motions will be so serious such that eq.(Q) may be far from the real situation. 
Moreover, there may be no coherent peak in the spectral function of one-particle Green 
function and we are unable to simply add a t'-term to the effective Hamiltonian. In other 
words, the t-term should play a more important role in this case. As has been shown in ref. 
|T0| , upon doping, the three- leg t — J ladder has two components - - a conducting Luttinger 
liquid coexisting with an insulating spin liquid phase. In order to discuss its low energy 
physics in the same spirit as the present paper, we need more understanding about the 
ground state of the undoped case and the behavior of holes in this ground state. Work along 
this line is under progress. 



APPENDIX A: BOSONIZATION RULES 



We list all bosonization rules we used in the following: 

^l{z) = 7= — exp{-i\/47i(f) L (z)}, 
\J lira 

^r{z) = -^=exp{iV4n(j) R {z)}, 
\/ lira 
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where z = r + ix, z is its complex conjugate and a is the short distance cut-off. iPl(z) and 
iPr{z) are left and right fermions, respectively. <Pl( z ) an d 4>r(z) are bosonic left and right 
movers, respectively. In terms of them, we define 0(r, x) = 4>l + 4>r and 6(r, x) = 4>l — 4>r. 
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